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Many extensions of the standard model predict the existence of hidden sectors that may contain 
unbroken abelian gauge groups. We argue that in the presence of quantum decoherence photons 
may convert into hidden photons on sufficiently long time scales and show that this effect is strongly 
constrained by CMB and supernova data. In particular, Planck-scale suppressed decoherence scales 
D oc to 2 /Mp\ (characteristic for non-critical string theories) are incompatible with the presence of 
even a single hidden U(l). The corresponding bounds on the decoherence scale are four orders 
of magnitude stronger than analogous bounds derived from solar and reactor neutrino data and 
complement other bounds derived from atmospheric neutrino data. 
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I. MOTIVATION 

Physics beyond the standard model typically predicts 
the existence of hidden sectors containing extra matter 
with new gauge interactions. There is no compelling 
reason why these extra sectors should be very massive 
if their interaction with standard model matter is suffi- 
ciently weak. This can be accomplished, e.g., in extra 
dimensional extensions with a geometric separation of 
sectors or in models where tree-level interactions are ab- 
sent and higher order corrections suppressed by the scale 
of messenger masses. 

The presence of these light hidden sectors can have 
various observable effects. For example, in the case of 
kinetic mixing [TJ[2] between a hidden U(l) and the elec- 
tromagnetic U(l), hidden sector matter can receive small 
fractional electro-magnetic charges. These mini-charged 
particles can have a strong influence on early Universe 
physics and astrophysical environments [3] • If the hidden 
sector U(l) is slightly broken by a Higgs or Stiickelberg 
mechanism we can have a situation analogous to neutrino 
systems with characteristic oscillation patterns between 
photons and hidden photons over sufficiently long base- 
lines [4]. 

Here we concentrate on another effect induced by the 
presence of these light hidden sectors: the sensitivity of 
photon propagation to sources of quantum decoherence, 
e.g. quantum gravity effects [5]. A heuristic picture de- 
scribes space-time at the Planck scale as a foamy struc- 
ture [H] , where virtual black holes pop in and out of exis- 
tence on a time scale allowed by Heisenberg's uncertainty 
principle [7J. This can lead to a loss of quantum infor- 
mation across their event horizons, providing an "envi- 
ronment" that might induce quantum decoherence of ap- 
parently isolated matter systems [5]. 

It is an open matter of debate whether quantum deco- 
herence induced by a quantum theory of gravity would 
simultaneously preserve Poincare invariance and local- 
ity [5J |7J |H1 GO HO] • A violation of energy and momentum 



conservation by particle reactions with a space-time foam 
could be reflected by an (energy dependent) effective re- 
fractive index in vacuo This could be tested, e.g., 
by the measurement of the arrival time of gamma rays 
or high energy neutrinos at different energies or by the 
propagation of ultra-high energy cosmic rays |12j . 

If, on the other hand, Poincare invariance is preserved, 
the presence of a non-trivial space-time vacuum can still 
be signaled by decoherence effects in systems of stable 
elementary particles [8 . A particularly interesting and 
well-studied case are neutrino systems, where the inter- 
play between mixing, mass oscillation and decoherence 
can influence atmospheric, solar and reactor neutrino 
data [T31 HU [TS], as well as flavour composition of as- 
trophysical high-energy neutrino fluxes |161 117) . 

If hidden sectors contain unbroken abelian gauge 
groups it is also feasible that the system of the elec- 
tromagnetic photon (70) and hidden photons ({74} with 
i > 1) experience energy and momentum conserving de- 
coherence effects: transitions between different photon 
"species" , 7 M — > 7„ , are allowed by gauge and Poincare 
invariance. This is the only alternative system to neu- 
trinos involving a stable and neutral elementary particle 
of the standard model, which can experience these de- 
coherence effects on extremely long (cosmological) time 
scales. 

The outline of this paper is as follows. We will start 
in Sect. |ll]with an outline of the Lindblad formalism of 
quantum decoherence. We will then discuss in Sect. |III| 
the effect of photon decoherence on the Planck spectrum 
of the cosmic microwave background and the luminosity 
distance of type la supernovae, respectively. This enables 
us to derive strong limits on various decoherence models. 



We comment in Sect. IV on the interplay of decoherence 
and photon interactions and outline a possible mecha- 
nism to extend the survival probability of extra-galactic 
TeV gamma rays. We finally summarize in Sect. [V] 
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II. QUANTUM DECOHERENCE 



III. OBSERVATIONAL CONSTRAINTS 



The Lindblad formalism is a general approach to quan- 
tum decoherence that does not require any detailed 
knowledge of the environment |18j . In the presence of 
decoherence the modified Liouville equation can be writ- 
ten in the form 



dp 

dt 



= -i[H,p]+V[p] 



(1) 



The Hamiltonian H can include possible background con- 
tributions, e.g. plasma effects for the photon. However, 
this does not effect the evolution of the density matrix in 
the absence of mixing between the gauge bosons. The de- 
coherence term D in the modified Liouville equation ([I]) 
can be written as 



V[P\ 
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[bj,pb l] + [bj p, 6 ] 



(2) 



where {bj} is a sequence of bounded operators acting 
on the Hilbert space of the open quantum system, 7i, 
and satisfying J2jb]bj £ B(H), where B(H) indicates 
the space of bounded operators acting on Tt. The dy- 
namical effects of spacetime on a microscopic system can 
then be interpreted as the existence of an arrow of time 
which in turn makes possible the connection with ther- 
modynamics via an entropy. The monotonic increase of 
the von Neumann entropy, S(p) — — Tr(plnp), implies 
the hermiticity of the Lindblad operators, bj = [19] . 
In addition, the conservation of energy and momentum 
can be enforced by taking [P M , bj] = 0. 

We will assume in the following that there is a total of 
N abelian gauge bosons 7 M with p, = 0, . . . , N — 1 includ- 
ing the photon j . The solution to Eq. (JTJ is outlined 
in Appendix [X] For simplicity, we assume degeneracy of 
the decoherence parameters Di — D which simplifies the 
photon survival probability after a distance x = t = L 
significantly (see App. |A|, 
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The energy behavior of D depends on the dimension- 
ality of the operators bj. We can estimate the energy 
dependence from gauge invariance. Possible combina- 
tions of the field strength tensors F^ u and G^ of the 
two U(l)s are [5] bj cx (F^G^Y cx (J . This restric- 
tion of the energy behavior to non-negative powers of 
uj may possibly be relaxed when the dissipative term is 
directly calculated in the most general space-time foam 
background [T7] . 

An interesting example is the case where the dissi- 
pative term is dominated by the dimension-4 operator 
b\ yielding the energy dependence D cx uj 2 /Mp\. This 
is characteristic of non-critical string theories where the 
space-time defects of the quantum gravitational "envi- 
ronment" are taken as recoiling D-branes, which generate 
a cellular structure in the space-time manifold |20j . 



In the following we will investigate the limits on quan- 
tum decoherence in the presence of hidden massless U(l)s 
from cosmological and astrophysical observations. Un- 
less otherwise stated, we will make the conservative as- 
sumption that there exists only a single hidden U(l) in 
addition to the standard model (N = 2) and that the 
decoherence rate of photons can be parametrized as 



D(z,u) = (l + z) p+1 D, 



LQP+ 1 

1\7p~ 

JW QD 



(4) 



where Mqp is the scale of quantum decoherence, not 
necessarily the Planck scale. We will consider the scale 
dependencies p = 0, ±1, ±2. 

The differential flux of photons from a source at red- 
shift z is reduced by the exponential factor 



. 1 1 

P 7 ^ 7 (z,w) = - + - exp 



d£'D(z',uj) 



(5) 



where the propagation distance I is given by 
dz = H(z)(l + z)d£ with Hubble parameter H . 
The Hubble parameter at redshift z is given by 

h 2 (z) = # 2 [(i - n m - n r ) + n m (i + zf + o.(i + z) 4 ] 

where Q, m h 2 ~ 0.128 and Q, r h 2 ~ 2.47 x 1CT 5 . The 
present Hubble expansion is Ho = h 100 kms -1 Mpc -1 
with ft ~ 0.73 EH- 



A. CMB Distortions 

In the standard big bang cosmology the cosmic mi- 
crowave background (CMB) forms at a redshift of about 
zcmb — 1100 after recombination of electrons and 
(mostly) protons in the expanding and cooling Universe. 
The CMB is well described by a Planck spectrum with a 
temperature of T = 2.725 ± 0.001 K [H], 



d 2 n 
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dudn 2ir 2 exp (f ) - 1 



(6) 



The high degree of accuracy (better than 1 in 10 4 around 
1 meV) between the CMB measurement and cosmological 
predictions is an ideal probe for exotic physics that could 
have effected the CMB photons in the redshift range < 
z < zcmb with energies from meV (z = 0) to eV (z = 
^cmb)- 

The influence of light particles coupling to the CMB 
has been studied previously for the case of axion-like par- 
ticles [53] , mini-charged particles [23] and massive hidden 
photons with kinetic mixing 25j. In the case of photon 
absorption or decoherence the observed spectrum is mod- 
ified as 
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FIG. 1: Modification of the CMB spectrum by decoherence 
effects. The dots show the residual CMB spectrum measured 
by COBE/FIRAS [22] (Uy (Janksy) = 1(T 26 W m" 1 Hz" 1 ). 
The upper panel shows the effect of D, = 0.01,0.1, 1 and the 
lower panel D, = 10' 4 . The curves show the accumulated 
deficit in photons compared to the "no decoherence" case. (A 
different shadowing is included to guide the eye.) 

As an illustration of the effect of decoherence, Fig. [T] 
shows the distortion of the CMB spectrum for the case 
p = 1 and A/qd = Mpi ~ 1.2 X 10 19 GeV for D* = 
0.01,0.1,1 (upper panel) and D* = 10 -4 (lower panel). 
The 3cr limits from a x^-fit of various models to the 
COBE/FIRAS data are shown in Tab. |T] (see also 
upper left panel in Fig. |3|. For comparison, the two 
right columns show limits on the decoherence parameter 
derived from solar and reactor neutrino data |15j and at- 
mospheric neutrino data (T3J [26]. We will comment on 
how to compare these bounds at the end of the section. 



B. SN Dimming 

Also the luminosity distance measurements of cos- 
mological standard candles like type la supernovae 
(SNe) [27j [28] are able to test feeble photon absorption 
and decoherence effects. The luminosity distance dz is 
defined as 

where C is the luminosity of the standard candle (as- 
sumed to be sufficiently well-known) and F the measured 



flux. In a homogeneous and isotropic universe this is pre- 
dicted to be 

*w -a +«)«.• (jf^). m 

with o,q 1 = H Q ^/\1- f2tot| and $k(f) = (sinh£,£,sin£) 
for spatial curvature k=— 1,0,1, respectively. If the 
photon flux of a source, located at distance z and ob- 
served in a (small) frequency band centered at w*, is at- 
tenuated by photon interactions or quantum decoherence 
the observed luminosity distance increases as 

d?\*) = dL{Z l ■ (10) 

The apparent extension of the luminosity distance by 
photon interactions and oscillations has been investigated 
in the context of axion-like particles [23 , 29 , hidden pho- 
tons [5U], chameleons [!JT] and mini-charged particles [3"2"] . 
One of the main attractions of these models is the possi- 
bility that the conclusions about the energy content of 
our Universe drawn from the Hubble diagram can be 
dramatically altered. We will briefly comment on this 
possibility at the end of this section. 

As an example, the upper panel of Fig. [2] shows the 
effect on quantum decoherence in the ACDM model as- 
suming a single hidden U(l) (N — 2). The luminosity 
distance of the SNe is shown as the difference between 
their measured apparent magnitude m and their known 
absolute magnitude M, 

m-M = 51og 10 d£ iM pc + 25. (11) 

As in the previous case we can derive 3ct limits for various 
decoherence models that are shown in Tab. [T] (see also 
upper right panel in Fig. [3j|. The limits on Planck-scale 
suppressed decoherence (p = 1,2) from SN dimming are 
slightly stronger than the CMB limits. 

However, we would like to emphasize that these limits 
depend on the cosmological model and the normalization 
of the SN data. If the SN dimming by photon decoher- 
ence is strong this can have an effect on the evaluation 
of cosmological data. We give an example in the lower 
panel of Fig. [2] where the decoherence effect could even 
be able to reproduce the observed SN luminosities from a 
CDM model. Note, that this model is not excluded by the 
corresponding CMB limits, though it is not compatible 
with the analogous bound from atmospheric neutrinos. 
In addition, one has to keep in mind that the photon 
frequency dependence of the parameter D causes a color 
excess (unless p = — 1) which is in conflict with obser- 
vation 28J. For the example shown in the right panel 
of Fig. [2] the color excess between the B and V band 
of the form, E[B - V] = A(m - M) B - A(m - M) v , 
is larger than 0.8 for redshifts 0.2 < z < 0.8, more 
than twice the median color excess observed from SNe at 
this distance [33]. Moreover, photon absorption as a SN 
dimming mechanism would violate the cosmic distance- 
duality, i.e. the luminosity and angular diameter distance 
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FIG. 2: Upper Panel: Enhanced SNe dimming by deco- 
herence effects assuming N = 2, p = 1, Mqd = Mp\ and 
£>, = 10- 6 ,10- 5 ,10 -4 assuming sources observed in a fre- 
quency interval centered at u t = 1 eV. The dots show the 
SNe la "union" compilation from Ref. |28j . The luminosity 
distance cLl is shown by the difference A(m — M) to an empty 
(fitot = 0) flat universe. The thin lines show the accumu- 
lated deficit in photons compared to the "no decoherence" 
case (standard ACDM model) indicated by the thick line. 
Lower Panel: As the upper panel, but now showing also a 
flat CDM model with Q m = 1 and Q A = (dotted line). For 
illustration, we consider a hidden photon model with p = 1, 
N = 2 and D» = 6 x 10 -6 and show the dimming for the B 
(A* ~ 440 nm, upper line) and V (A* ~ 550 nm, lower line) 
band. Whereas the overall dimming effect is practically in- 
distinguishable from the ACDM model, the strong reddening 
of the starlight from the energy dependence of D oc tv 2 /Mpi 
is challenged by the data |33| . 



relation d^^A = (1 + z) 2 (34]. Hence, it is unlikely that 
decoherence can fully account for SN dimming. Never- 
theless, it could have an effect on the evaluation of cos- 
mological data. 

The dependence on Mqd of the upper bound on for 
various models of D as derived from CMB distortions and 
SN dimming are summarized in the top panels of Fig. [3] 
The lower panels show analogous limits derived from so- 
lar and reactor neutrino data |15j as well as atmospheric 
neutrino results [T31 [25]. It is apparent from the plots 
(as the intersection of the bounds ) that these probes are 
based on different energy scales, meV-eV in the case of 
the CMB, eV in the case of SNe, multi-MeV in the case of 
solar and reactor neutrino data and multi-TeV in the case 



of atmospheric neutrinos. 1 Consequently we see that for 
models with p < the bounds from from CMB distor- 
tions and SN dimming are the strongest. One must how- 
ever bear in mind that while bounds from neutrino data 
are "unavoidable" once the quantum decoherence effects 
exist, the effects on photon propagation discussed here 
rely on the extra assumption of the presence of massless 
hidden [/(l)'s. Besides this, we also note that, generi- 
cally decoherence effects that grow with energy are bet- 
ter constrained from higher energy atmospheric neutrino 
data. For example for the case D oc uj 2 /Mqd (p = 1 in 
our notation) the possible decoherence of photons in the 
presence of a single hidden U(l) is about four (five) or- 
ders of magnitude stronger constrained from CMB (SNe) 
data than the corresponding limit from solar and reac- 
tor neutrino data |15j (c/. T ab. |T| but the bound derived 
from atmospheric neutrinos [26 is much stronger. 

The dimensionless parameter appearing in our def- 
inition of the decoherence scale Q is naturally expected 
to be of order 1. Hence, atmospheric neutrino data al- 
ready severely constraints Planck-scale suppressed deco- 
herence of the form D = D^oj 2 /Mp\. Similarly, our anal- 
ogous limits on photon decoherence in the presence of 
hidden U(l)s require <C 1 for p = 1. In reverse, a 
quantum theory of gravity predicting Planck-scale sup- 
pressed decoherence is not compatible with the existence 
of unbroken hidden U(l)s unless, unnaturally, <C 1. 

Higher order quantum decoherence with p > 2 is so 
far only weakly constrained by neutrino systems or by 
photon/hidden photon systems in the meV to eV energy 
range. We will speculate in the following section, how 
high energy gamma ray sources could possibly explore 
this unconstrained parameter space. 

IV. PHOTON PROPAGATION 

Decoherence effects can also have interesting effects on 
the propagation of photons in the presence of photon 
absorption. We can account for photon absorption effects 
in the Liouville equation ([TJ by a contribution 

V r [p] = -l{U ,p}, (12) 

where T = 1/A is the photon absorption rate, 2 e.g. in the 
intergalactic photon background (BG) via 7 + 7bg — > 
e + + e~, and Ho is the photon projection operator. 

The photon survival probability in the presence of a 
single hidden U(l) can be readily solved from the modi- 



1 The sensitivity to astrophysical neutrinos in the TeV-PeV region 
have also been discussed in Ref. [ 16). This probe could signifi- 
cantly improve the bounds on decoherence for the cases p = 1, 2. 

2 Since we consider point-source fluxes in the following we will treat 
the reaction 7 + 7bg ~~ * e+ + e ~ as an absorption process of the 
photon. Subsequent electro-magnetic interactions of secondary 
er^ will contribute to the diffuse GeV-TeV background. 
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atmospheric v 


p = -2 


9.5 x 10 


3.9 x 10 


8.3 x 10 


1.0 x 10 


D — — 1 


9.7 x 10~ 66 


3.4 x 10~ 62 


8.8 x 10~ 44 a 




p = 


3.9 x 10~ 35 


2.9 x 10" 34 


1.0 x 10" 22 


1.2 x 10" 27 c 


p = \ 


1.3 x 10~ 5 


2.3 x 10 -6 


0.11 a 


1.6 x 10~ 12 c 


p = 2 


3.5 x 10 23 


1.8 x 10 22 




910 c 



a [l5] 3ct C.L. from their Fig. 1 with n = p + 1 and 70/GeV = D*(M QD /GeV) p . 
i> |13| The authors consider the case p = — 2 with /j 2 = 2D, A/ 2 D 

c 26 Bounds are at the 90% C.L. Their notation corresponds to n = p + 1 and = D*(M p i/AfQD) p . 



TABLE I: The 3cr limits on D* derived from the COBE/FIRAS data [22] and SN data [28] for a decoherence scale given by 
Eq. Q. We assume N = 2 and Mqd = Mpi. Note that the limits are in general stronger for iV > 2 by up to a factor 2. For 
comparison the last two columns show the limits from reactor and solar [T3] , and atmospheric neutrinos [131 126] . 
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FIG. 3: 3cr C.L. upper bound on D* as a function of Mqd for various models of D. We show the bounds on from the CMB [22] 
(upper left panel) and from type la SNe [28 (upper right panel) for various models p — 0, ±1, ±2 of the decoherence parameter 
D = D,lo p+1 /Mq D . Darker colors indicate regions which are excluded by multiple models. Note that the limits become 
generally stronger for W > 2 by up to a factor 2. For comparison, the lower plots show the limits from solar and reactor 
neutrino data [T^] (3a C.L.) and atmospheric neutrino data, [TS] (3cr C.L.) and [5B] (90% C.L.). 
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fied Liouville equation and gives (see App. |A|for details) 



P 7 ^ 7 = e -^ {D+T) 



cosh ( ^VD 2 +T 2 



, r sinh ( - Vd 2 + r 2 



(13) 



We can estimate the sensitivity of extra-galactic TeV 
gamma ray sources to decoherence effects as 



(14) 



In particular for Mqd ~ Mpi , ~ 10 for p = 1 or 
-D* ~ 10 -6 for p = 2. Hence, this probe has the potential 
to be more sensitive to Planck-scale suppressed decoher- 
ence than existing neutrino data 1 ,'! 2(i by many orders 
of magnitude. We will discuss in the following possible 
signals of decoherence in the spectra of TeV gamma ray 
sources. 

Figure [3] shows the survival probability for three dif- 
ferent values of D. The functional behaviour can be 
easily understood as follows: The hidden U(l) serves as 
an invisible "storage" of photons during propagation. If 
DA > 1 decoherence quickly equalizes the number of 
photons and hidden photons. This results into a decrease 
of the photon survival probability to 1/2 (for N = 2) for 
L < A. On the other hand, inelastic scattering only af- 
fects photons. Therefore at larger L photons can be "re- 
plenished" by the decoherence of the unabsorbed hidden 
photons into photons. Hence, for propagation distances 
L > A the presence of hidden photons increases the pho- 
ton survival probability. 

For the general case of N U(l)s and strong decoherence 
D>rwe can express the photon survival probability at 
a distance much larger than the decoherence scale L 3> 
l/D as (see App. |A| for details) 



P 

7->7 



1 

N 



exp 



NX J 



(15) 



This can have an important effect on the spectra of TeV 
gamma ray sources. 3 Firstly, if the onset of decoherence 
appears at energies covered by the spectra one could ob- 
serve a step-like drop of the flux by a factor XJN. And 
secondly, as long as the absolute source emissivity of pho- 
tons is unknown (such that the pre-factor 1/N in Eq. ( 15 ) 



3 Note that in the absence of hidden sectors, as long as ui < v (v 
being the scale of electroweak symmetry breaking ) the photon 
survival probability is 1 independently on whether the Higgs po- 
tential around virtual black holes has its minimum at the trivial 
vacuum v = 246 GeV or at the unbroken vacuum v = 0. For 
photon energies beyond the electroweak breaking scale, one may 
theorize over possible decoherence effect between the neutral SM 
gauge bosons, tn this case the conservation of energy and mo- 
mentum in decoherence effects, which forbids transitions of the 
form 7 — * Z, is a crucial assumption. 
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FIG. 4: Photon survival probability in the presence of deco- 
herence with a single hidden U(l) (N = 2). The black line 
shows the photon attenuation, e.g. in the inter-galactic pho- 
ton background via 7 + 7bg — > e + + e~, in the absence of 
decoherence (D — 0). The dashed curve shows decoherence 
with DX = 1 and the dotted curve DA 2> 1. 



gets renormalized) , the expected spectral cut-offs of the 
sources could be shifted according to an extend photon 
interaction length A c ff — NX. These effects can be clearly 
seen in Fig. [4] for the case N = 2. For a TeV gamma ray 
source at 100 Mpc and Planck-scale suppressed decoher- 
ence this requires > 10~ 27 (D» > 10 -11 ) for p = 1 
(p = 2). In comparing this with Tab. [I] and Fig. [3] we see 
that there is ample room for models that could have such 
an effect on the spectra. 



V. CONCLUSIONS 

We have discussed the effects of quantum decoherence 
of photons in the presence of hidden sector U(l)s. Quan- 
tum decoherence in the system of photon and hidden pho- 
tons could be induced by a foam-like structure of space- 
time in a quantum theory of gravity, where virtual black 
holes pop in and out of existence at scales allowed by 
Heisenberg's uncertainty principle. 

We have shown that these decoherence effects are 
strongly constrained by the absence of photon disappear- 
ance in the cosmic microwave background. Furthermore 
quantum decoherence as an additional source of starlight 
dimming can also be constrained by the luminosity dis- 
tance measurements of type la supernovae. Consequently 
based on the standard ACDM model we can derive con- 
straints on the decoherence in the presence of hidden 
U(l)s. In principle, the decoherence effect could be 
strong enough to influence the evaluation of cosmolog- 
ical data. However, color dependencies in the dimming 
via photon decoherence are not favored by the data and 
could be used to derive further constraints. 

Our main results are summarized in the upper panel of 
Fig. [3] and compared with the present bounds from non 
observation of decoherence effects in neutrino systems. 
For models with decoherence parameter D constant or 
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decreasing with energy, the bounds from from CMB dis- 
tortions and SN dimming are the strongest. One must 
however bear in mind that while bounds from neutrino 
data are "unavoidable" once the quantum decoherence 
effects exist, the effects on photon propagation discussed 
here rely on the extra assumption of the presence of 
massless hidden U(l)s. We also note that, generically 
decoherence effects which grow with energy are better 
constrained from the higher energy atmospheric neutrino 
data. 

We have also discussed the interplay between photon 
absorption and decoherence. This effect can become im- 
portant for distant TeV gamma ray point-sources if the 
decoherence length is much smaller than the photon in- 
teraction length. Assuming N — 1 additional U(l)s this 
could leave characteristic features in gamma ray spectra 
in the form of step-like drops by factors l/N or by an ef- 
fective increase of the absorption length A to A ff = NX. 
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APPENDIX A: LINDBLAD FORMALISM 

We outline the solution to the Liouville Eq. in the 
presence of quantum decoherence. The density matrix 
p and Lindblad operators bj can be expanded in a basis 
of hermitian matrices that satisfy the orthonormality 
condition Tr (F^F V ) — <5 M „/2. Without loss of generality 



we consider a basis with (-Fo)ij 
we have 



P 



Explicitly, 



(Al) 



cone" coordinates x = (x—t) /2 and t = (x+t)/2. In these 
new coordinates Eq. fan can be written dp/di = T>[p\. 

Hence, the coefficients of the free equations of motion 
satisfy the differential equation 



d_ 

Of 



/V = - X D »»P>> > 



with D^o = A)p 



and 



hi ■ 



(A2) 



(A3) 



fe,/, 



where fijk are structure constants defined by [Fi,Fj] — 

* ^2k fijkFk- 

The solution of d t po — is trivial and requires po(t) = 
pa = 2Tr(n Q F ) = ^2/N for all species a 
diagonalizable by a matrix M, (M~ 1 DM) i j - 
can write the final solution as 



If Dij is 

DiSij, we 



p — 



1 

N 



(A4) 



This reduces to Eq. (J3j> in the case Di = D using 2 = 
p 2 ^ following from p 2 (0) — p(0) and Tr[p] = 1. 

We can extend the Liovillc Eq. (JT|) by a photon decay 
term of the form (12 1. The modified equation can be 



readily solved in the case of N = 2, taking Fi — \<Ji with 
Pauli matrices <Ji. In this basis the photon projection 
operator in the term ( 12 1 has the form n = F +F 3 . And 



the general solution of the photon survival probability is 
given in Eq. (13 1. 



In the case of strong decoherence, i.e. at propagation 
distances L ^> 1/D and D ^> T, we can derive the 
asymptotic solution (15 1 of the general survival proba- 



bility with N U(l)s in the following way. In the presence 
of strong decoherence we can assume that at any step 
during the evolution Trfn^p] ~ Trfn^p] and, in particu- 
lar, Tr[p] ~ iVTrplop]. Since the trace of the decoherence 
term vanishes we can derive the asymptotic differential 
equation 



Tr[p] ~ NTr[U Q p] ~ -rTr[n p] 



(A5) 



The free propagation of photons and hidden photons 
(H = —id/dx) can be readily solved in terms of "light- 



This has the solution (15 1, noting that P 7 _^ 7 = Trpop] 



and initially Tr[n o p(0~ l/N 
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